REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite 
1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. _ 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  TYPE  REPORT  3.  DATES  COVERED  (From  -  To) 

14-09-2007  Conference  Paper  15  June  07 


4.  TITLE  AND  SUBTITLE  5a.  CONTRACT  NUMBER 


Statistical  characteristics  of  microwave  signals  scattered  from  a  randomly  rough 
surface 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

AFRL/SNHE 
80  Scott  Drive 

Hanscom  AFB  MA  01731-2909 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 

61102F 


5d.  PROJECT  NUMBER 

2304 


5e.  TASK  NUMBER 

HE 


5f.  WORK  UNIT  NUMBER 

01 


8.  PERFORMING  ORGANIZATION 
REPORT 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Electromagnetics  Technology  Division  SC:  437490 

Sensors  Directorate 
Air  Force  Research  Laboratory 
80  Scott  Drive 

Hanscom  AFB  MA  01731-2909 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

DISTRIBUTION  A.  Approved  for  public  release;  distribution  unlimited. 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

AFRL-SN-HS 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

AFRL-SN-HS-TP-2007-0003 


13.  SUPPLEMENTARY  NOTES 

Supported  by  Air  Force  Office  of  Scientific  Research  (AFOSR)  for  AFRL/SNHE,  IN-HOUSE  work  unit  2304HE01 . 
Published  in  Proceedings  of  IEEE  International  Symposium  on  Antennas  and  Propagation,  Honolulu  HI,  June  2007. 
Cleared  for  Public  Release  by  ESC/PA:  ESC  06-1478  dtd  19  Dec  06.  This  is  a  work  of  United  States  Government  and  is 
not  subject  to  copyright  protection  in  the  United  States. 


14.  ABSTRACT 

This  paper  is  an  investigation  of  the  statistical  characteristics  of  the  amplitude  of  signals  scattered  from  a  randomly  rough 
surface.  When  the  area  of  illumination  and  the  grazing  angles  are  large  the  amplitude  of  scattered  signals  are  found  to 
obey  Rayleigh  statistics.  On  the  other  hand,  when  the  grazing  angles  are  very  small  and  the  resolution  size  is  small  it  has 
been  found  Lognormal  statistics  form  a  good  fit.  However,  in  the  intermediate  domain  we  find  that  no  particular  distribution 
has  overall  best  fit.  The  conclusion  is  based  on  a  statistical  analysis  of  the  simulated  data  generated  by  rigorous  rough 
surface  scattering  analysis  along  with  Monte  Carlo  simulation  of  sample  surfaces.  The  scattering  calculations  are  based 
on  integral  equation  formulation  of  scattering  and  are  solved  using  method  of  moments. 


15.  SUBJECT  TERMS 

statistics  of  amplitude  of  scattered  signals,  randomly  rough  surface,  Monte  Carlo  simulation,  non-Rayleigh  statistics, 
rigorous  wave  scattering  computations 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

Unclassified 


b.  ABSTRACT 

Unclassified 


c.  THIS  PAGE 

Unclassified 


17. LIMITATION 

18. NUMBER  OF 

OF  ABSTRACT 

PAGES 

uu 

6 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Saba  Mudaliar 


19b.  TELEPHONE  NUMBER  (include  area  code) 

N/A 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


11 


Statistical  Characteristics  of  Microwave  Signals  Scattered  from  a 
Randomly  Rough  Surface 

Saba  Mudaliar1  and  Freeman  Lin2 

1  Air  Force  Research  Laboratory,  Hanscom  AFB,  MA  01731-2909,  USA 
E-mail:  saba.mudaliar@hanscom.af.mil 
2  ARCON  Corporation,  Waltham,  MA  02451,  USA 


Introduction 


Scattering  of  microwave  signals  from  rough  surfaces  has  been  extensively  studied  [1-3]. 
However,  a  vast  majority  of  such  studies  is  on  the  characteristics  of  the  normalized  scatter¬ 
ing  cross  section.  This  paper  is  an  investigation  of  the  statistical  characteristics  of  scattered 
signals.  When  the  area  of  illumination  is  large  it  is  seen  that  the  scattering  from  a  rough 
surface  may  be  considered  to  be  contributions  from  a  large  number  of  individual  cells.  One 
may  hence  appeal  to  the  central  limit  theorem  and  show  that  the  amplitude  of  the  scattered 
signal  satisfies  the  Rayleigh  distribution.  This  observation  was  first  made  by  Rayleigh  a 
hundred  years  ago  for  the  general  case  of  randomly  scattered  signals.  Beckmann  demon¬ 
strated  that  similar  arguments  and  deductions  may  be  employed  for  the  case  of  scattering 
from  rough  surfaces.  It  is  apparent  that  certain  conditions  arc  necessary  to  deduce  that 
the  amplitude  statistics  are  Rayleigh  distributed.  When  such  conditions  do  not  exist  the 
amplitude  statistics  arc  not  Rayleigh  and  in  such  situations  we  would  like  to  know  which 
statistical  distribution  is  most  representative.  It  has  been  found  that  the  K-distribution  [4]  is 
a  good  model  for  randomly  scattered  signals.  We  consider  this  and  some  other  appropriate 
non-Rayleigh  statistical  distributions.  For  our  study  we  use  numerically  simulated  signals 
scattered  from  a  randomly  rough  surface.  We  employ  the  moment  method  to  estimate  the 
parameters  of  the  various  distributions.  We  next  use  Kolmogorov-Smirnov  statistic  to  de¬ 
termine  which  of  the  distributions  most  closely  fit  the  simulated  data. 

Description  of  the  Problem 


We  have  a  perfectly  conducting  randomly  rough  surface  given  as  z  =  k(x)  which  is  planar 
on  the  average.  (  is  a  zero-mean  stationary  random  process  independent  of  y.  For  TE  case 
the  incident  field  is  taken  as  E,  =  yE-„  where 


Ej  =  exp  \ik(x  sin  0,  -  z  cos  #;)[!  +  w(r)]}  exp 


w 


r(x  +  ztan9,)u 


(1) 


w  denotes  the  width  of  the  Gaussian  beam,  and 


(dr)  = 


1 

(kw  cos  6j)2 


-^(x  +  ztanF,)2  -  ll. 
w-  j 


(2) 


Equation  (1)  describes  a  Gaussian  beam  often  used  in  numerical  simulations.  For  the  TM 
case  we  have  H,  =  yH,  where  //,  is  identical  to  the  £,  in  (1).  Our  main  interest  is  in  studying 
the  statistical  characteristics  of  the  fields  scattered  by  the  rough  surface. 
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Formulation 


The  reduced  wave  equations  for  time-harmonic  waves  arc 

V2E  +  k2  E  =  0  V2H  +  k2  H  =  0,  (3) 

where  k  is  the  free  space  wave  number.  Since  the  surface  is  independent  of  y  the  polarization 
of  the  incident  waves  is  preserved  on  scattering.  Therefore  E  =  yE  and  H  =  yH.  Thus  we 
have  the  following  scalar  Helmholtz  equations  for  TE  and  TM  cases  respectively, 

V2£  +  k2E  =  0  V2H  +  k2H  =  0.  (4) 

The  boundary  condition  on  the  PEC  surface  is  z  x  E  =  0  on  the  surface  S .  This  implies  that 

E(rs)  =  0  dnH(  r,)  =  0  rseS,  (5) 

for  TE  and  TM  cases,  where  dnH  is  the  normal  derivative  on  the  surface  S .  Employing 
the  Green’s  theorem  we  convert  the  differential  equation  system  into  the  following  integral 
equation  system: 

E(r)  =  Ej(r)  -  f  dr'G(r',r)dn,E(r')  H(r)  =  H,(r)  +  f  dr'H(r')dn,G(r',r).  (6) 

Js  Js 

The  essential  task  is  to  evaluate  the  surface  fields  dnE  and  H  which  arc  governed  by  the 
following  integral  equations: 


dnE(rs)  =  2dnE,{rs)  -2  £  dr'd„G(r',  rs)dn'E(r'),  (7) 

H(rs)  =  2 Hi(rs)  +  £  dr'H(r')dn'G(r',rs).  (8) 

The  bar  across  the  integral  sign  denotes  that  the  integral  is  a  principal  value  integral.  We 
employ  the  forward-backward  algorithm  [5,6]  to  obtain  numerical  solutions  to  the  above 
integral  equations.  Inserting  these  solutions  in  (6)  we  calculate  the  bistatic  scattering  coeffi¬ 
cients.  Proceeding  thus  we  carried  out  numerical  simulations  based  on  Pierson-Moskovitz 
surface  spectrum  for  various  surface  conditions,  beam  widths,  and  angles  of  incidence  and 
observation. 


Statistical  Characteristics  of  Scattered  Waves 

Simulated  results  thus  obtained  arc  valuable  for  understanding  the  nature  and  characteris¬ 
tics  of  various  asymptotic  approximations.  Further,  the  procedure  provides  comprehensive 
statistics  of  scattered  signals.  One  important  quantity  of  interest  is  the  probability  density 
function  (pdf)  of  the  amplitude  of  scattered  signals.  The  pdf  of  signals  scattered  from  a  ran¬ 
dom  collection  of  scatterers  has  been  investigated  by  Rayleigh.  For  the  case  of  a  randomly 
rough  surface  this  problem  was  studied  by  Beckmann.  The  idea  is  to  formulate  the  problem 
as  a  random  walk  process  and  assume  that  the  surface  is  composed  of  contributions  from 
a  large  number  of  independent  individual  cells.  From  this  one  can  deduce  that  the  distri¬ 
bution  of  amplitudes  of  scattered  signals  is  Rayleigh  distributed.  However,  this  derivation 
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is  based  on  several  assumptions  and  when  they  arc  violated  we  depart  from  Rayleigh  sta¬ 
tistics.  There  are  several  non-Rayleigh  statistical  distributions  that  have  been  postulated 
for  such  situations.  Some  of  these  arc  lognormal,  Weibull,  K-distribution,  and  modified 
K-distributions.  The  modified  K-distributions  arc  extensions  of  the  K-model,  obtained  by 
using  different  functions  to  describe  the  texture.  For  our  study  we  have  considered  the 
gamma  texture,  log-normal  texture  and  the  inverse  gamma  texture.  The  pdfs  of  the  statisti¬ 
cal  models  cited  above  arc: 

Log-normal: 


m  = 


:V2 


:  exp 


Tier 


2a2 


ln  I  - 


for  s  >  0 


Weibull: 


K-distribution: 


m  = 


cs 


c — 1 


exp 


for  s  >  0 


f(s)  = 


V^v/f 

2v_1T(v) 
Generalized  K  -  gamma  texture: 


for  5  >  0 


2b  r  ( 
f(s)  ~  Uv)  ( 


v\vh  r 
M/  Jo 


r.Vb-2 


expi--(-Tl  \dT 


Generalized  K  -  Log-normal  texture: 


m  = 


V2 


r>o 

=  f 

2na2  Jo 


exp 


—  -  T^-ln 2[^-\\dT 
t  2a-  \2m) 


Generalized  K  -  Inverse  gamma  texture: 


m  = 


2  s/3T{a  +  1) 


(J3s2  +  l)^+1T(a) 


The  Rayleigh  distribution  has  a  one  parameter  pdf.  All  the  rest  have  two  parameter  pdfs. 
We  use  the  first  two  moments  to  determine  these  parameters  from  the  simulated  data  using 
Pearson  moments  method.  Let  9\  and  62  be  the  two  parameters  and  let  pi  and  /Ji  be  the  first 
two  moments  for  a  particular  model.  If  sn  is  the  amplitude  of  the  simulated  scattered  signal 
corresponding  to  the  n-th  realization,  then  the  moments  may  be  calculated  as 


1  Nr 

-Y 


mk  =  T7 

Nr 


where  Nr  is  the  number  of  realizations  considered  for  the  study.  Then 

Hk(d\,02)  -  mk  k  =  1,2 

is  the  pair  of  equations  to  be  solved  to  determine  the  parameters. 


4811 


This  is  a  fairly  simple  and  straightforward  procedure.  However,  there  is  no  information  as 
to  whether  the  parameters  thus  obtained  will  be  unbiased  and  efficient.  Fisher’s  maximum 
likelihood  estimation  method  guarantees  efficiency  although  it  is  not  as  simple  to  implement 
as  the  moments  method.  Let  the  Nr  variate  density  function  be  fyr.  Then  the  likelihood 
function  is  fyr (xj,  .V2,  •  •  •  ,  .v.vjd] ,  6b)-  The  parameters  6\  and  (h_  are  determined  by  solving 
dgl  fNr  -  0  and  dg2f^r  =  0.  The  next  obvious  question  is  which  pdf  most  closely  fits  our 
simulated  data.  To  compare  the  data  with  a  particular-  model  it  is  more  convenient  to  use 
the  cumulative  distribution  function  (cdf).  Denoting  the  proposed  cdf  as  F(s),  and  that 
composed  by  using  the  computed  data  as  Fc(s),  we  calculate  the  Kolmogorov-Smirnov 
statistic  D  =  sup  \F(:(s)  -  F(s)\.  To  gauge  the  pdf  under  consideration  we  calculate  the 
probability  that  D  is  larger  than  the  observed  value,  which  is  given  as 

Prob  (D  >  observed  value)  =  Q  ||  y[N\-  +  0.12  +  — —  j  D 


where 

OO 

Q(t)  =  2  ^(-ly-V2^. 

7=1 

Calculating  this  probability  for  the  various  proposed  distributions  enables  us  to  choose  the 
one  that  best  fits  the  data  for  the  particular  situation.  On  carrying  out  this  procedure  we 
find  that  the  pdf  that  fits  the  data  best  varies  with  angles  of  incidence  and  observation, 
beam-width  of  illumination  and  surface  characteristics.  Only  for  angles  close  to  normal  is 
the  Rayleigh  pdf  appropriate.  There  is  no  one  pdf  that  suits  all  situations.  We  find  that 
different  regions  and  situations  require  different  pdfs.  Certain  pdfs  are  more  suitable  for 
certain  situations  than  others. 
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